Theory of Mott insulator/band insulator heterostructure 
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A theory of heterostructures comprised of LaTiOs (a Mott insulator) and SrTi03 (a band in- 
sulator) is presented. The band structure of the Ti d-electrons is treated with a nearest neighbor 
tight-binding approximation; the electric fields arising from the La 3+ /Sr 2+ charge difference and 
the carriers are treated within a Hartree approximation; and the on-site interactions are treated 
by unrestricted Hartree-Fock. The phase diagram as a function of interaction strength and layer 
number is determined and predictions are made for optical conductivity experiments. A note worthy 
finding is that the edges of the heterostructure are generally metallic. 



I. INTRODUCTION 

"Correlated electron systems" (such as transition 
metal oxides) are materials in which strong electron- 
electron or electron-lattice interactions produces be- 
havior incompatible with the standard "density func- 
tional plus Migdal-Eliashberg" theory which describes 
most compounds. The past decade has seen tremen- 
dous progress in the physics and materials science of 
correlated-electron systems. Significant improvements in 
crystal and film growth, in measurement techniques and 
in theory have led to a much improved understanding 
of the bulk properties of these materials. An impor- 
tant finding is that correlated electron systems exhibit a 
multiplicity of interesting phases (superconducting, mag- 
netic, charge and orbitally ordered) along with associated 
novel excitations. For recent reviews, see Ref. 01 or the 
articles in Ref. 0- 

The recent success in treating bulk properties suggests 
that the time is ripe for a systematic study of the sur- 
face and interface properties of correlated electron sys- 
tems. In addition to its basic importance as a fundamen- 
tal question in materials science, correlated electron sur- 
face/interface science should provide the necessary scien- 
tific background for study of potential devices exploiting 
correlated electron properties, because essentially any de- 
vice must be coupled to the rest of the world via motion 
of electrons through an interface, and for study of cor- 
related electron nanostructures, because essentially the 
defining property of a nanostructure is a high surface 
to volume ratio. The fundamental interest of bulk cor- 
related electron materials lies in the novel phases they 
exhibit, and we therefore suggest that the fundamental 
issue for the nascent field of "correlated electron surface 
science" is " how does the electronic phase at the surface 
differ from that in the bulk" ; in other words, "what is 
the electronic surface reconstruction." 

This question has begun to attract experimental at- 
tention. Hesper and co-workers have shown that the 
[111] surface of K 3 Ceo differs from bulk because of charge 
transfer caused by a polar surfaced Matzdorf and collab- 
orators have demonstrated that in the correlated elec- 
tron system Cao.gSro.iRu03 (which exhibits Mott metal- 
insulator transition), the surface layers remain metal- 



lic down to a lower temperature than does the bulk 
system. 4 Izumi and co-workers have fabricated "digital 
heterostructures" composed of different transition metal 
oxides and have demonstrated changes in electronic phase 
and other properties depending on the thicknesses of dif- 
ferent layers^ In an experimental tour-de-force, Ohtomo, 
Mullcr, Grazul, and Hwang have demonstrated the fab- 
rication of atomically precise digital heterostructures in- 
volving a controllable number n of planes of LaTi03 (a 
correlated-electron Mott-insulating material) separated 
by a controllable number m of planes of SrTi03 (a more 
conventional band-insulating material) and have mea- 
sured both the variation of electron density transverse 
to the planes and the dc transport properties of the 
heterostructure^ Their work opens the door to controlled 
studies both of correlated electron physics in confined di- 
mensions and of the behavior of the interface between a 
correlated system and a uncorrelated one. 

Many physics and material science issues arise in con- 
sidering the behavior of correlated electrons near surfaces 
and interfaces. Atomic reconstruction may occur, and 
may change the underlying electronic structure. For ex- 
ample, the authors of Ref. y| argue that a change in tilt 
angle of the surface RuOg octahedra increases the elec- 
tronic hopping, thereby allowing the metallic phase to 
persist to lower T. Also, as noted e.g. by Hesper and 
collaborators, a change in structure will lead to changes 
in Madelung potentials, and to the screening which helps 
define the values of many-body interaction parameters^ 
"Leakage" of charge across an interface may change den- 
sities away from the commensurate values required for in- 
sulating behavior. Substrate induced strain is well known 
to change the behavior of filmsii 

Sorting out the different contributions and assessing 
their effect on the many-body physics is a formidable 
task, which will require a sustained experimental and 
theoretical effort. The experiment of Ohtomo et al. of- 
fers an attractive starting point. In this system, the near 
lattice match (1.5 % difference in lattice parameter) and 
chemical similarity of the two components (LaTiOa and 
SrTiOa) suggests that atomic reconstructions, strain, and 
renormalizations of many-body parameters are of lesser 
importance, so the physical effects of electronic recon- 
struction can be isolated and studied. Further, the near 
Fermi surface states are derived mainly from the Ti d- 



orbitals^ and correspond to narrow bands well described 
by tight-binding theory. However, the orbital degeneracy 
characteristic of Ti provides an interesting set of possible 
ordered phases. 

In this paper we undertake a theoretical analysis of 
the correlated electron behavior to be expected in lattice- 
matched digital heterostructures of the type created by 
Ohtomo et al& We focus on electrons in the Ti-derived 
bands, and include the effects of the long-ranged electric 
fields arising both from the La atoms and the electronic 
charge distribution. We treat the extra on-site interac- 
tions via a Hartree-Fock interaction. We calculate the 
electronic phase diagram as a function of on-site inter- 
action parameter and number of La layers and for the 
relevant phases determine the spatial variation of charge, 
spin and orbital densities. We obtain a complex set of 
behaviors depending on interaction strength and number 
of La layers. Generally, we find a crossover length of ap- 
proximately three unit cells, so that units of six or more 
LaTi03 layers have a central region which exhibits bulk- 
like behavior. The outmost ~ 3 layers on each side are 
however metallic (in contrast to the insulating behavior 
of bulk LaTiOs). For very thin superlattices the ordering 
patterns differ from bulk. We calculate optical conduc- 
tivity spectra and show that this is a revealing probe of 
the electronic structure. 

The rest of this paper is organized as follows: Section 
II defines the model, parameter values, and method of 
study. Section III presents our results for the T = 
phase diagram as a function of interaction strength and 
number of LaTiC>3 layers. Section IV presents our results 
for the spatial dependence of the charge density in the 
small U regime where neither spin nor orbital ordering 
occur and provides, for this simpler case, an overview of 
the general features of the electronic structure. Section 
V studies the onset of spin and orbital order as the in- 
teraction and layer thickness are increased. Section VI 
discusses in more detail the special case 1 layer. Section 
VII study the interrelation between the metallic behavior 
and the subband structure, and a stability under parame- 
ter variation. Section VIII presents representative results 
for the optical conductivity, and shows how these may be 
used to elucidate the electronic structure. Finally, sec- 
tion IX presents a summary of our findings and important 
future directions, and implications for the general ques- 
tions of correlated electron surface and interface science, 
and is designed to be useful to readers uninterested in 
the details presented in sections II- VIII. 
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FIG. 1: Schematic figure of the model used in the present 
study. Open and crossed circles show the positions of Sr and 
La ions. (La layer number n — 2) The dots show the positions 
of Ti ions. x,y axes are chosen to be parallel to the La plane, 
and z axis is perpendicular to the plane. 



aiding in the growth of high quality digital heterostruc- 
tures. 

In this paper we consider an infinite crystal of SrTiC>3, 
into which n adjacent [001] planes of LaTiC>3 have been 
inserted perpendicular to one of the Ti-Ti bond direc- 
tions, as shown in Fig. ^ We choose the z direction to 
be perpendicular to the LaTiC>3 planes, so the system has 
a (discrete) translation symmetry in the xy direction. 

The relevant electronic orbitals are derived from 
the Ti i2g-symmetry d-states, and maybe labeled as 
d X y, d xz , d yz . The derived bandsAS are to a good approx- 
imation given by a nearest neighbor tight binding model 
with hopping parameter of magnitude t ~ 0.3 eV and 
spatial structure given by the Slater-Koster formula, 11 
so the d xy states disperse only in the xy plane etc. 

We take the electrons to hop between Ti sites accord- 
ing to a nearest neighbor hopping and to feel a potential 
defined by (i) the Coulomb force arising from the extra 
charge on the La relative to Sr (ii) the Coulomb force 
arising from the electrons on other Ti sites and (Hi) on- 
site "Hubbard U and J" interactions with other electrons 
on the same site. We take the form of the on-site interac- 
tions determined by Mizokawa et alJ^ and adopt values 
as discussed below. Thus 
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II. FORMALISM 

Both LaTi03 and SrTi03 crystallize in the simple 
ABO3 perovskite structure^ (more precisely, very small 
distortions occur which we neglect here) and as noted by 
Ref. the lattice constants of the two materials are al- 
most identical; aL a Tio 3 — 0-SrTiO 3 — 3.9 A, minimizing 
structural discontinuities at the interface and presumably 



and similarly for xz,yz. We have 

#cL = V c {Ri)p d (Ri) 
with p d (Ri) = £„ CT d\ ia d aia and 
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Here the Ri are the positions of the Ti( B ) sites in the 
ABO3 lattice and R^ a label the actual positions of the 
La ions, which reside on the A sites. 

We denote the dielectric function of the host lattice 
by e. An interesting issue arises here: SrTi03 is a nearly 
ferroelectric material^ The static dielectric constant be- 
comes very large at long wavelength and low temper- 
atures, but e is much smaller at high frequencies, room 
temperature, or short length scales. Also the polarization 
P will increase more slowly at higher fields, and relevant 
quantity is P/E. In this paper we have chosen e — 15 as 
a compromise between these effects. We discuss the con- 
sequences of different choices of e in the conclusion. We 
emphasize that incorporating the ferroelectric tendencies 
of SrTi03 (including the associated lattice distortions) in 
a more realistic manner is a important question for future 
research. 

Finally, the onsite H is 

H cm-sit e = U^n ai] n ali + {U' - J) ^ naiaUMa 
a a>b,<j 

a^b a^b 

+ J 'Y. d U d w d lii dbi i- ( 5 ) 

a/b 

The critical issue is the strength of the on-site repul- 
sion. For definiteness we follow other studies which em- 
ploy the ratios V — 7U /9 and J = U/9 which are very 
close to those determined by Mizokawa^ Many work- 
ers have used the value U ~ 5-6 eV ~ 18i-20i estimated 
from high energy spectroscopies. 8 However, optical con- 
ductivity studies of LaTiC>3 and related compounds such 
as YTiC>3 find rather small gaps, in the range 0.2-1 eV^ 
suggesting U ~ 2.5 eV ~ 8t. In view of this uncertainty 
we investigate a range U from ~ 6t-20t. 

To study the properties of Htot, Eq. l(T[l. we employ 
the Hartree-Fock approximation replacing terms such as 
by n aia (n bia ) + (n alrT )n bi a] orbitally off-diagonal 
expectation values {{d ai(T dbia)) of the type considered by 
Mizokawa and Mochizuki are stable only in the presence 
of a GdFeC>3 type distortion which we do not consider. 
While not exact, the approximation reveals the correct 
trends and in particular reveals insulating, ordered states 
in the parameter regimes where these exist. To imple- 
ment the Hartree-Fock approximation we assume an ini- 
tial distribution of site, spin and orbital occupancies, 
obtain the one electron potential by factorizing the in- 
teraction terms as described below, compute the band 
structure, obtain new densities, and iterate until conver- 
gence is obtained for some [/-values. Many iterations 
(~ 10 3 ) are required to obtain well-converged solutions, 
essentially because of the delicate balance of particle dis- 
tributions needed to screen the long range part of the 
Coulomb interaction. 

There are two types of solutions to the one-electron 
equations: bound states, which decay as \z\ — » 00, and 



continuum states, which do not. As usual in heterostruc- 
ture problems the bound states give rise to sub-bands, 
some of which are partially occupied. The ground state is 
obtained by filling the lowest sub-bands up to the appro- 
priate chemical potential (determined by charge neutral- 
ity); the interaction-related terms in H@ oul and H^l-site 
are then recomputed and the procedure is repeated until 
self-consistency is obtained. Charge neutrality requires 
that the total density of electrons in the bound-state sub- 
bands equals the total density of La ions. However, the 
interplay between electron-La attraction and electron- 
electron repulsion leads (in almost all of the cases we 
have studied) to a very weak binding of the highest-lying 
electron states; indeed for large U the Fermi level of the 
partially filled sub-bands is only infinitcsimally below the 
bottom of the continuum bands. Also in some regions of 
the phase diagram, there are several locally stable solu- 
tions and it is necessary to compare energies, to deter- 
mine the ground state. Between 200 and 3000 iterations 
were required, with the larger number needed either when 
the Fermi level adjoins the continuum states or when sev- 
eral different states are very close in energy. 

We now mention some specific details of our solution to 
the Hartree-Fock equations. We seek solutions which are 
bound in the z directions; i.e. functions xy) which 
decay as layer index i — > ±00, and which depend on the 
transverse (xy) coordinates as discussed below. 

For each orbital a, the functions f a (i] xy) obey a single 
particle Schodinger equation which we write as 

t«a 2 D 2z [ip a (v,xy)) + H^<p a (i;xy) = E<p a (i;xy). (6) 

Here D^zf is the discretized second order derivative 

tpi+i ~ 2ipi + <fii-i and t" = t = 0.3 eV for a = xz, yz 

(i) 

and = for a = xy. The term Hxy contains the dis- 
persion in the plane of the layers, as well as terms aris- 
ing from the long-ranges Coulomb interaction [Eq. 
and terms arising from the Hartree-Fock approximation. 
These terms depend on the on-site interactions and on 
the charge, spin and orbital density in layer i. A gen- 
eral mean-field state would break translation invariance 
in the xy plane, in a manner whose amplitude would de- 
pend on layer index i, thereby mixing the motion in x, y 
and z directions, leading to a very numerically involved 
three dimensional self-consistency problem. In order to 
keep the computations within reasonable bounds and ex- 
plore wide ranges of parameters, we have in most calcu- 
lations restricted attention to state which preserve trans- 
lational invariance in the xy plane. For these states the 
xy and z motions decouple, and physical quantities may 
be computed from the solution of a Id problem combined 
with an integral over a two-dimensional momentum. We 
also performed some investigations of states with a two- 
sublattice xy plane structure. In this case, after Fourier 
transformation on the in-plane coordinate the in-plane 
Hamiltonian may be written 

H (i) = (Vi+A* e p \ 



with Vi the part of the potential which is independent of 
sub-lattice, A' the part which alternates between sublat- 
tices and e p the in-plane dispersion in the reduced Bril- 
louin zone. The variation of the mixing term A 1 with 
layer requires a separate treatment, lengthening the nu- 
merical analysis considerably. 



III. PHASE DIAGRAM 

In this section we present and discuss the calculated 
phase diagram. The geometry of the heterostructure dis- 
tinguishes between the xy orbital and the {xz, yz} or- 
bitals. As noted above, we assume translation invariancc 
in the xy plane. The symmetries which may be broken 
are therefore spin rotation and orbital rotation (xz and 
yz) and inversion along with translation invariance in the 
xy plane. Also, although it does not have precise mean- 
ing except in the n — ► oo limit, it is physically sensible to 
interpret results in terms of ^-direction translation sym- 
metry breaking when the spin or orbital density oscillates 
with respect to the total charge density. 

Our calculated phase diagram is shown in Fig. |3 For 
reasons of computational convenience in scanning a wide 
range of parameters, we considered mainly phases with 
translation invariance in the xy plane, however for n = oo 
and n = 1 we considered also an xy-plane two sublattice 
symmetry breaking. We found in agreement with previ- 
ous calculation a 12 ! 15 ! 16 ! 17 that the fully staggered phase is 
favored at n — oo, but xy-plane symmetry broken states 
could not stabilized in the one layer case. We have not 
yet studied more general symmetry breakings for inter- 
mediate layer numbers, but physical arguments presented 
below strongly suggests that these phases only occur for 
larger numbers of layers (n > 6) . 

Four phases are found: a small U phase with no bro- 
ken symmetry, and intermediate U phase with in-plane 
translation-invariance spin order, but no orbital order, 
and a large U phase with both spin and orbital order. 
The lower U transition line varies smoothly with layer 
number, while the larger U transition is essentially inde- 
pendent of layer number for n > 1. Further, the n = 1 
intermediate U phase is ferromagnetic whereas for n > 1 
the intermediate U phase is antiferromagnetic. As will 
be shown in sections below this behavior can be related 
to the sub-band structure but the essential reason is that 
for small n the charge density is spread in the z direction, 
so no layer has a density near 1. 

For all n, the lower-f/ transition is found to be second 
order. As U is increased above this value, the magnetiza- 
tion increases rapidly and the upper transition is strongly 
first order. Fig. [3] shows an example of the procedure 
used to determine the location of the phase boundary 
and the order of the transition. The upper panel of fig- 
urc|s3displays the energies of several different phases. We 
identify U c \ as the point where the different energies con- 
verge (this is most easily done from a plot, not shown, of 
the energy differences on a expanded scale) and U C 2 from 
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FIG. 2: Theoretical ground state phase diagram as a function 
of the on-site Coulomb interaction U and inverse of the La 
layer number n. Parameter values are chosen to be U' = 
7(7/9 and J = 17/9. e 2 /(eat) — 0.8 which corresponds to 
e ~ 15 with a lattice constant a = 3.9 Aand transfer intensity 
t — 0.3 eV. The triangle is the critical U to the (ir, tt, 7r) 
orbital ordering at n = oo case. The broken line shows the 
exDected Dhase transition to the (n.n.n) antiferromasnetic 
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FIG. 3: (a) Energy as a function of U (n — 1). Open squares: 
the spin and orbital-disordered state, filled squares: the spin- 
order orbital-disorder state, open circles: the spin and orbital- 
ordered state, (b) Square of magnetization as a function of U 
in the intermediate coupling region (n = 1). 



the lowest crossing point of the energy curves. That the 
energies of two qualitatively different states cross at U c i 
identifies this as a first order phase transition, that the 
energies merge at U c \ suggests a second order transition. 
Further evidence is provided by panel (b), which shows, 
on an expanded scale, the magnetization as a function of 
U — Ud, and suggests a continuous decrease to zero. 

The comparison to the n — oo limit is subtle. In bulk 
LaTi03 exhibits a (71", tt, tt) type antiferromagnetic or- 
dering. Theoretical calculations (apparently confirmed 
by very recent NMR experiment, and x-ray and neu- 
tron diffraction experiments)i2*iS suggest a four sublat- 
tice structure which is very close to a (0, 0, 7r)-type orbital 
ordering 20 differing slightly from the (0,0, tt) ordering 
studied here. Stabilizing the observed state apparently 
requires a lattice distortion not included in the model 
studied here. As U is increased from zero the n — > 00 
limit of the model considered here has a phase transition 
which we believe to be of second order to an incommen- 
surate antiferromagnetic state with a wave vector which 
is an extremal spanning vector of the Fermi surfaces of 
the bands arising from two of the orbitals (say xz, yz) 
and which turns out to be very close to (0, 0, 7r). (In fact, 
for reasons of numerical simplicity we studied (0, 0, tt) 
ordering and found a very weakly first order transition.) 
This transition is followed by a strongly first order tran- 
sition to one of a degenerate manifold of states charac- 
terized by ferromagnetic spin order and (tt, tt, tt) orbital 
order. To maintain continuity with the heterostructure 
calculations we have also suppressed the (tt, tt, tt) order- 
ing and located the phase boundary to the (metastable) 
(0, 0, 7r)-orbital spin-ferromagnetic state. We believe the 
(tt, tt, 7r)-orbital spin-ferro state we have found is a rea- 
sonable surrogate for the actual Mochizuki-Imada state 
found in experiment. Although this state is favored at 
71 — > 00, the physical origin, explained below in more 
detail, of our inability to stabilize states with broken in- 
plane translation invariance at n = 1 suggests that the 
(tt, tt, tt) phases only occur for rather thick superlattices. 
The essential point is that, for n < 6, the solution in the 
large U limit consists of several partially filled sub-bands, 
which have effectively minimized their interaction energy 
but which gain considerable kinetic energy from motion 
in the xy plane. Breaking of xy translation symmetry 
would reduce this kinetic energy gain without much de- 
creasing the already saturated interaction energy while 
z-direction kinetic energy is quenched by the confining 
potential. 

The estimate of U = 8t-10t derived from optical mea- 
surements on bulk, and from the recent work on the na- 
ture of the ordered phase suggests an interesting series of 
phase transitions may occur as layer thickness is varied. 



IV. SMALL U REGION 

In this section we present and discuss the results ob- 
tained for U < U c i where no symmetry breaking is 
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FIG. 4: Calculated total electron density in the small U region 
(U = 6t). Our results (not shown) for larger U values are 
almost identical. 



observed. Although the system geometry implies xy 
and {xz,yz} orbitals are inequivalent, the occupancies 
of those orbitals are found to be almost the same; the 
ratio at the center cite is given by (n xz — n xy )/n xz = 
0.023,0.098 and 0.049 for n = 1,2 and 3, respectively, 
at U — 6t. Therefore, we mainly focus on the spatial 
distribution of the total electron density below. 

Numerical results for the electron density for several 
choices of La layer n and U = 6t are shown in Fig. 0] 
Here we set the origin of z axis at the center of n La lay- 
ers. Electron density is symmetrically distributed around 
z = 0. The density at the center site rapidly increases 
with increasing n, and exceeds 0.9 at n = 6. However, for 
this and larger U values, the charge density never exceeds 
unity, even at the center sites. Fig. 0] shows clearly that 
for thick layers the length scale over which density vary 
is about three unit cells, and this number is almost in- 
dependent of n, and depends somewhat on e and weakly 
on [7i2i Thus, at least 6 layers seem to be required be- 
fore bulk behavior is recovered in the central region. As 
is well known, SrTiOa is close to ferroelectric, and the 
long wave-length, and linear response dielectric constant 
becomes about 20000 at the lowest temperature^ The 
present results involve e at short scales, and are clearly 
outside the linear response regime. However, tempera- 
ture dependent changes in the density profile should oc- 
cur and are worth experimental study. Another inter- 
esting feature to note is that for all e and U studied the 
electron and hole distributions at large n are almost sym- 
metric at the edge z — n/2; n to t(z) and 1 — n tot (n/2 — z) 
fall on almost the same line (not shown). 

It might be interesting to compare the present nu- 
merical results (n — 1,2) with the experimentally ob- 
served Ti 3+ spatial distributions^ The observed Ti 3+ de- 
cay length of Ti 3+ is about 2 nm (5 unit cells) in both 
n = 1 and 2, which is longer than the present theoret- 
ical results (about 3 unit cells). The reported La-ion 
distribution suggests that the spatial resolution of the 



EELS measurement is about 0.4 nm (1 unit cell). Inclu- 
sion of such effect makes the density distribution profiles 
closer to the experimental ones, but the experimental dis- 
tributions are still broader than the present theoretical 
ones. One possibility is La/Sr intcrdiffusion (not con- 
sidered here). An alternate origin could be a chemical 
shift of the Ti t2 g level. Here we have considered only 
the long-range Coulomb energy, but an additional con- 
tribution arising from the slightly different local structure 
[small ionic radius of La 3+ (1.045 A) than that of Sr 4+ 
(1.13 A)}22 could bring the La-Ti d-levels closer to Sr-Ti 
d-levels. This effect has not been included in this study. 



V. ONSET OF POLARIZATION; LARGER U 

This section presents results for large U values, where 
spin and orbital orderings are observed. The main focus 
is on how the orderings are developed. 

Figure[S](a) and (b) show the electron densities per or- 
bital in the intermediate coupling regime (U = St) with 
n = 2 where only spin ordering is observed. The differ- 
ence between the {xz, yz} and xy occupancy arise only 
from the symmetry. As a comparison, electron densities 
per orbital and spin in the metastable paramagnetic state 
are also shown (crossed squares). From the difference 
between the up and down electron occupancy, one can 
notice that the spin density oscillates around the center 
site, (although no site is fully polarized) In contrast for 
n = 1 and 8t < U < 12t, the magnetic state is different: 
the spins align ferromagnetically. The behavior at n = 1 
can be attributed to the lower charge density at each 
site implying low occupancy of each sub-band, so a small 
shift between the up and down electron sub-bands can 
stabilize the fully polarized state. On the contrary, fer- 
rimagnetic spin ordering at n > 2 comes from a delicate 
balance between the confinement potential and on-site 
Coulomb repulsion. 

The electron densities per spin and orbital in the strong 
coupling regime (U = lOt) are shown in Fig. [S] (c) and 
(d) . The difference between xz and yz orbital occupancy 
clearly indicates antiferro-type orbital ordering along the 
z direction (0, 0, 7r). One can also notice that spin polar- 
ization is almost complete. A tiny amount (less than 1 %) 
of down spin electrons exist only in the xy orbitals. Note 
that the orbital ordering in the strong coupling regime 
is different from that expected in the bulk limit where 
translational symmetry in xy plane would also be broken. 
In the small n region, the symmetry in xy is expected to 
be conserved because this symmetry is favorable to gain 
the kinetic energy by two dimensional electron motion. 



VI. n = 1, POSSIBILITY OF STAGGERED 
IN-PLANE ORDERING 

In this section, we discuss special features of the n = 1 
case. Within the Hartree-Fock approximation we have 




FIG. 5: Electron density per spin and orbital in the interme- 
diate coupling regime U = St (a,b) and the strong coupling 
regime U = Wt (c,d) for a heterostructure with two La layers 
(n — 2). Open and filled symbols show the densities of the 
up and down spin electrons, respectively. Circles and squares 
in (a,c) are for the xz and yz, respectively. Crossed squares 
in (a,b) show the electron density per spin and orbital in a 
metastable paramagnetic state. 



not found any stable states which break xy plane trans- 
lational invariance. We see from Fig. 0] that at essen- 
tially all U the occupancy of any particular orbital state 
is low. For U > 13t the orbital ordering in the z direc- 
tion (say xy occupancy to the left of the La layer and yz 
occupancy to the right), along with the ferromagnetism 
effectively eliminates the on-site interaction contribution 
to the Hartree-Fock energy, while the fact that the one 
electron is shared between the two bands means that at 
each sub-band there is a large kinetic energy gain of or- 
der — — sin pp per sub-band. We believe this physics ex- 
plains why we at n — 1 have been unable to to stabilize 
a state with two-sublattice in-plane symmetry breaking. 
In bulk at large U, the kinetic energy is quenched and the 
(it, 7T, 7t) state is stabilized relative to the (0, 0, it) state by 
a superexchange of approximately 2t 2 / (V — J) ~ 3t 2 /U 
but in the n = 1 superlattice the band energy gain is ap- 
proximately — ^ sin pp. Indeed this argument suggests 
that until the band filling becomes close to 1 (which only 
happens for n ^ 6) the kinetic energy gain in rry-invariant 
states is larger than the superexchage energy gained by 
ordering so that for a range of small n (perhaps n < 6) 
staggered in-plane ordering may be neglected. 

We do expect the xy-translation invariant states we 
have found to be unstable to some form of weak in- 
commensurate density wave ordering, because the re- 
sulting sub-band structure has a one-dimensional char- 
acter from the xz and yz sub-bands. However, this or- 
dering is expected to be weak because the residual in- 
teractions, not already included in our Hartree-Fock so- 
lution, are not strong. For example, in the U > 13t 
regime, the on-site interactions are fully quenched by 



the ferromagnetic ordering so incommensurate charge 
density wave is driven only by the 2p F component of 
the long ranged Coulomb interaction, which we esti- 
mate from the sub-band charge distribution with half 
width £ to be V c (2p F ) m C (2e 2 /ea) with C = 
/ °° duexp[-u 2 (7r£/a) 2 ]/Vl + u 2 w 0.3-0.4. For the pa- 
rameters used here V c X2p F ~ Xln(Ep /T) with A <C 0.1, 
so the instability occurs at a very low T. We note, how- 
ever, that the precise critical behavior of the second order 
small U transition found in our Hartree-Fock approxima- 
tion is likely to be complicated by the additional effects 
of 2pp instabilities in the quasi one dimensional bands. 
We leave this interesting issue for future research. 



VII. METALLIC EDGE AND STABILITY 
UNDER PARAMETER VARIATION 

The experimental results of Ohtomo et al. 6 show that 
the interface between the two types of insulators sup- 
ports a metallic state. From numerical results for thick 
layer (Fig. |3J, one finds a central region where n ~ 1, in 
which one expects bulk-like behavior, and an edge region 
with n <C 1, which one might expect to exhibit metallic 
behavior. In this section, we study the interrelation be- 
tween the subband structure and metallic behavior. We 
show that the edge region is responsible for the metal- 
lic behavior, and study the dependence of the metallic 
behavior on parameter changes. 

Within the Hartree-Fock method we use here, physics 
of the metallic edge is manifested as follows. There 
are many bound-state solutions (wave functions decay- 
ing as \z\ is increased away form the heterostructure), 
whose dispersion in the in-plane direction gives rise to 
sub-bands. For thin heterostructures, all sub-bands are 
partially filled (implying metallic behavior in the het- 
erostructure plane) whereas for thick heterostructures 
(in ordered phases), some sub-bands are fully filled and 
some are partly filled. The fully filled sub-bands have 
z-direction wave functions implying charge density con- 
centrated in the middle of the heterostructure, whereas 
the partially filled bands have charge concentrated near 
the edges. The upper panel of Fig. shows the den- 
sity profiles and the occupancy of partially filled bands 
for the heterostructure with n = 6 and U/t = 10 in a 
orbitally ordered ferromagnetic phase. We observe par- 
tially filled bands corresponding to metallic behavior at 
the heterostructure edges where the density drops from 
- 1 to - 0. 

Two key assumptions underlying the calculations pre- 
sented so far are the translation invariance of the on- 
site energy Eq (because in both components of the het- 
erostructure the electronically active ion is octahedrally 
coordinated Ti) and relatively large value for the dielec- 
tric constant (to represent the nearly ferroelectric nature 
of SrTiOz). We investigate the sensitivity of this result 
to changes in these parameters. The lower panel of Fig.0 
shows the changes in the density profile and in the oc- 




FIG. 6: Upper panel: density profile (dashed line) and den- 
sity in partially occupied bands of heterostructure with n = 6, 
U/t = 10 and e = 15 calculated in orbitally-ordered ferromag- 
netic phase. Lower panel: Density profile (dashed lines) and 
density in unoccupied bands (light and heavy solid lines) for 
on-site binding to Ti sites in heterostructure (on site energy 
for Ti) between two La planes eo; on site energy for "edge" 
Ti sites (between Sr and La planes) eo/2. 



cupancy of partially filled bands occurring if the on-site 
energy is changed for sites near La sites. Here, we have 
chosen the site energy of those Ti sites which sit between 
two Sr planes to be zero, the energy of a Ti site be- 
tween two La planes to be Eq < and the energy of a 
Ti site between a La plane and a Sr plane to be Eq/2. 
For £o = several partially filled bands exist, leading to 
a high density of metallic electrons. As £o is decreased, 
the electronic structure rapidly rearranges so that most 
bands become fully filled. However, even for the largest 
£q the geometry ensures that one electron always remains 
to be shared between two nearly degenerate bands, one 
on either side of the heterostructure, so the metallic be- 
havior remains robust. 

We may similarly consider changes occurring as the di- 
electric constant is varied. Fig.[7|shows the paramagnetic 
phase density profile occurring as the dielectric constant 
is varied over a wide range. The main effect is to decrease 
the tailing of charge density far into the SrTiO^, layer. 
However, there again remains one electron to be shared 
between the two edge subbands, the basic phenomenon 
of partly filled bands at the edge is unaffected. 

Summarizing, for reasonable parameters we find that 



s 




FIG. 7: Dependence of total charge density on dielectric con- 
stant calculated for paramagnetic phase of 6 layer heterostruc- 
ture with U/t = 10 in paramagnetic phase. 



there are always partially filled bands, corresponding to 
metallic behavior at the heterostructure edges. The pos- 
sibility of metallic behavior at the edge of a Mott insu- 
lator system appears to be confirmed by the experimen- 
tal results of Ohtomo et al. 6 would be interesting study 
further both experimentally and by more sophisticated 
theoretical methods. 
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VIII. OPTICAL CONDUCTIVITY 

There are two classes of solutions to the Hartree- 
Fock equations: bound states, which have wave functions 
which decay as \z\ is increased away from the La layers, 
and continuum states, extended in all three directions. 
Although the existence of sharply defined single-particle 
states at all energies is an artifact of the Hartree-Fock 
approximation, we expect this qualitative structure to 
survive in a more sophisticated treatment. Therefore, in 
this section we use the Hartree-Fock approximation to 
show how optical conductivity can provide information 
on the nature and filling of the bound states. We also 
present a qualitative discussion of the effect of interac- 
tions beyond Hartree-Fock. 

We consider electric field directed perpendicular to the 
La-planes (parallel to z) and use the Peierls-phase ap- 
proximation to determine the current operator 

Jz = -it [4i<Ai-M<r - h.c] . (8) 

ia—xz.yz 

We then evaluate the usual Kubo formula using the 
Hartree-Fock eigenstates. Note that in the nearest neigh- 
bor tight binding approximation the d xy states do not 
couple to z-electric field, so we will not discuss them in 
this section. 

Figure. [S] shows the calculated T = conductivity for 
the case of one La layer, and two values of U: U = Qt 
(light curve), for which the system is in the spin and 
orbital disordered phase, and U = lOt (heavy curve), 



FIG. 8: (a) Optical conductivity for n = 1 with U = 6t (small 
U regime) and lOt (intermediate U ferromagnetic state). Elec- 
tric field E is chosen to be parallel to z. Delta functions 
are approximated as Lorentzian with the half-maximum half- 
width O.lt. (b) Schematic illustrations of sub-band energy 
structure at U — 6t (left) and U = lOt (right). Shaded ar- 
eas show the continuum state. Broken lines show the Fermi 
energy. 



for which the system is in the ferromagnetic, orbitally 
disordered state. 

For U = 6t the xz, yz part of the spectrum (shown 
in the lower left panel) consists of one four-fold degener- 
ate bound-state sub-band, containing approximately 1/6 
electron per spin per orbital, as well as empty continuum 
bands. The optical spectrum therefore consists of one 
feature, corresponding to a bound-continuum state tran- 
sition. The absorption is peaked at the bound state en- 
ergy, the width reflects the degree of overlap between con- 
tinuum and bound state energies. The oscillator strength 
is related to the kinetic energy in the z direction, and is 
approximately 10 4 17 _1 cm _1 /La. We expect this spec- 
trum to be very little affected by correlation effects be- 
yond Hartree-Fock, because the basic bound state ener- 
getics are fixed by charge neutrality and Coulomb inter- 
action, and the final states are continuum states, which 
are delocalized in space. 

Also shown in the upper panel of Fig. [S] is the spec- 
trum corresponding to n = 1 and U — lOt. For these 
parameters the electronic structure has changed (lower 
right panel of Fig. 0; there are now two sub-bands; one 
more strongly bound than in the U = 6t case and holding 
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more electrons (but still only approximately 1/3 filled) 
and one only very slightly filled band. The calculated 
o{ui) correspondingly exhibits three features, a weak low 
energy (u) ~ O.lt) feature arising from transitions from 
the very slightly filled bound to the continuum, a sharp 
higher energy (ui ~ 0.9£) feature arising from the inter- 
subband transition (allowed because the lower sub-band 
is even under z «-> — z and the upper is odd) and a very 
broad feature from lower bound-continuum transitions. 
The intersubband transition contains about 2/3 of the 
total oscillator strength; the remainder is mainly in the 
bound-state to continuum peak. We note that as ex- 
pected on general grounds, the presence of a bound-state 
to bound-state transition sharply reduces the weight in 
the continuum (compare U = St, Wt), essentially because 
of the requirement that the continuum states be orthog- 
onal to both bound states. 

Within the model we have used the intersubband tran- 
sition is a delta function, because lifetime effects have 
been neglected and the x, y and z direction dispersions 
decouple, so the 2 sub-bands disperse in the same way as 
k x or k y is varied. The dispersions decouple because we 
have adopted a nearest neighbor tight binding model, 
however this is believed to be reasonably accurate in 
practice. A potentially more significant source of broad- 
ening is the "Hubbard {/" interaction. Experience from 
simpler models suggests that for these U -values (of order 
of the critical value for the bulk system Mott transition) 
and electron concentration (less than 0.4 electrons per 
layer, far from the Mott value n = 1) these effects are 
not too severe: essentially because the probability of 2 
particle collisions scales as n 2 , only a small fraction of 
the spectral weight will be shifted to an "upper Hubbard 
band" feature at an energy ~ U , and lifetime broadening 
will be rather less than the Fermi energy, so roughly we 
expect the peak to remain unchanged. These arguments 
apply also to the other cases we consider below. 

We now discuss the n = 2 case. In the small U param- 
agnetic, orbital disordered case, the electronic spectrum 
exhibits two bound-state sub-bands; arising from a more 
strongly bound even parity state and a weakly bound odd 
parity state. The band structure is very similar to that 
shown in the right panel of Fig. 0(b), but without spin 
polarization. The optical conductivity (not shown) is es- 
sentially the same as the U = lOt curve shown in Fig. [S] 
Increasing U to U = 8t leads to a spin ordered (here ferri- 
magnetic) orbitally disordered state, with four sub-bands 
[three below fi and partially occupied] and correspond- 
ingly two bound-state to bound-state transitions (one for 
each spin) and three bound-state to continuum transi- 
tions (the middle one of which is evidently extremely 
weak) . 

As a final example, we show the optical spectrum and 
energetics for the spin-ferro orbital-antiferro U = Wt 
state. Five bound states now occur, and correspondingly 
three sharp peaks (recall the optical selection rule al- 
lows coupling only from even to odd parity states) . Four 
bound-state-continuum transition should be visible (the 
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FIG. 9: Optical conductivity for n = 2 with U — St (interme- 
diate U ferrimagnetic state) and lOt (large U ferromagnetic 
orbital ordered state). Schematic illustrations of sub-band en- 
ergy structure at U — 8t (b) and U = Wt (c). Shaded areas 
show the continuum state. Broken lines are the Fermi energy. 
Remark: band curvature not shown. 



highest-lying xz even parity state is essentially unoccu- 
pied) but again all are weak and the lowest-lying (arising 
from the odd d xz band) is especially weak because the 
low-lying even-parity continuum states must be orthogo- 
nal to the even-parity bound state, which absorbs all the 
oscillator strength. 

To summarize: the electronic spectrum consists of a 
number of different sub-bands (with the detailed struc- 
ture depending on number of La layers and on the na- 
ture of the ordering, if any). The optical spectrum for 
z-polarized light is predicted to consist of relatively sharp 
peaks, which may be assigned to the different optically 
allowed intersubband transitions, along with broad fea- 
tures (often rather weak) associated with transitions from 
sub-bands to the continuum. 



IX. SUMMARY AND CONCLUSION 

In this paper we studied theoretically the phase di- 
agram and electronic properties of a "correlated het- 
erostructure" model involving n layers of a material 
which in bulk is a Mott insulator, embedded in an in- 
finite band insulator. The specific features of the model 
we study were chosen to reproduce the LaTiOa/SrTi03 
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superlattice system studied by Ohtomo et ah, but we 
hope that our results will shed light also on the more 
general question of the physics of the interface between 
a strongly correlated and weakly correlated systems. 

A crucial feature of the experimental LaTiOa/SrTiOs 
system studied by Ohtomo et al. is the almost perfect 
lattice match between the two systems. These authors 
argued that this implies that the only difference between 
the Mott insulating and band insulating regions arises 
from the different charge of the La(3+) and Sr(2+); in 
particular the crystal structure and atomic positions are 
expected to remain relatively constant throughout the 
heterostructure. Of course, the asymmetry present at 
the LaTiOs/SrTiOs interface must induce some changes 
in atomic positions: a Ti06 octahedron is negatively 
charged, and so if it sits between a Sr plane and a La 
plane it will be attracted to the latter, and also distorted, 
because the positively charged Ti will tend to move in 
the opposite direction. The experimentally determined 
Ti-Ti distances shown in Fig. 1 of Ref. |g, along with the 
distortion in that paper, suggests that the changes in Ti- 
Ti distance are negligible. In this circumstance, changes 
in O position along the Ti-O-Ti bond change hoppings 
only in second order. We therefore neglected these effects 
and assumed that the electronic hoppings and interaction 
parameters remain invariant across the heterostructure. 
However, we emphasize that properly accounting for the 
effect of atomic rearrangements inevitably present at sur- 
face and interface is crucial. We further note that lattice 
distortions appear to be important in stabilizing the ob- 
served bulk state, but may be pinned in a heterostruc- 
ture. Extending our results to include these effects is an 
important open problem. 

In the calculation reported here the heterostructure 
is defined only by the difference (+3 vs +2) of the La 
and Sr charge. The calculated electronic charge density 
is found to be controlled mainly by electrostatic effects 
(ionic potentials screened by the electronic charge distri- 
bution). Results, such as shown in Fig. 0] for U = 6t, 
are representative of results obtained for a wide range of 
on-site interaction U. We find generally that significant 
leakage of charge into the band insulator region occurs. 
The width of the transition region must depend on the 
relative strength of the z-direction hoppings and the con- 
fining potential. For the parameters studied here, the 
transition region is about 3 layers (so one needs about 7 
La layers to obtain a central region with bulk behavior) , 
and that the n — 1 case is special because the La counter- 
ions sit between the Ti planes. We note, however, that 
the moderate to large [/-values which we study do en- 
sure that even for thick superlattices the charge density 
on the central layers never becomes greater than unity. 

The calculated width is found to be somewhat less than 
that found by Ohtomo et al., broadening our distribution 
by the experimental uncertainty (obtained from the La 
positions) leads to calculated width about 2/3 of mea- 



sured ones. Whether the difference arises because we 
have overestimated the confining potential or for some 
other reason such as La/Sr intcrdiffusion remains to be 
determined. 

We now turn to our calculated phase diagram shown in 
Fig. |21 It is expected on general grounds that decreasing 
the number of La layers will raise the interaction val- 
ues required for obtaining ordered states. Further, the 
specific structure of the system of present interest im- 
plies an odd-even alternation. Both of these features are 
indeed observed in our calculation. Our calculation, in 
combination with the U -values inferred from optics, sug- 
gests that the La/SrTi03 system may be in the inter- 
esting {U ~ 7t-12t) parameter regime in which one or 
more phase boundaries may be crossed by varying La 
layer number. A further point, perhaps especially rel- 
evant to systems such as titanates which allow for or- 
bital ordering and orbitally dependent hopping, is that 
the structural anisotropy intrinsic to the heterostructure 
may favor different ordering patterns than those found 
in bulk. Thus we find for the thinnest heterostructures 
different phases with more translational invariance than 
found in bulk. We note, however, that in the Ti case 
the "edge states" have a quasi one-dimensional character 
with an incommensurate filling, perhaps favoring incom- 
mensurate charge, spin or orbital ordering at the edge. 

The issue of the transport properties of the het- 
erostructure is an important open question, especially 
in light of the interesting transport data of Ref. 0. A 
crucial experimental finding is that metallic conduction 
is always observed, consistent with our prediction of con- 
ducting edge states. However, we do not have a qual- 
itative understanding of the Hall data, because within 
the Hartee-Fock approximation the Hall resistivity is ob- 
tained by adding up the contributions of the different 
sub-bands. The contributions to o~ xy arising from the 
quasi one-dimensional xz, yz sub-bands are controlled by 
weak deviations from the nearest neighbor hopping ap- 
proximations, which control the mixing between the two 
bands and the reconstruction of the Fermi surface of the 
crossing point. It is difficult to make general statements 
without more detailed band structure information, but 
some degree of compensation is expected. 

Other important issues for future research include 
study of correlation effects beyond the Hartree-Fock ap- 
proximation, inclusion of the coupling between atomic 
rearrangements, orbital ordering and electronic hopping 
parameters, and more sophisticated treatment of the di- 
electric properties of SrTiOs. 
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